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I.  INTRODUCTION 


The  deactivation  of  vibrationally  excited  HF  by  the  products  of  die 
F + H2  reaction  directly  affects  the  performance  of  the  HF  chemical  laser. 

The  deactivation  of  HF(v)  by  H has  been  studied  by  both  theoretical  calcula- 
tions  ' and  experimental  investigations.  " The  trajectory  calculations  of 
Thompson*  and  Wilkins2  were  the  first  indications  that  H atoms  might  be  ef- 
ficient deactivators  of  the  upper  levels  of  HF.  Subsequently,  in  two  experi- 
3 6 

mental  studies,  ’ rate  coefficients  for  HF(v  = 3)  deactivation  were  found 

13  3 

that  were  greater  than  10  cm  /mol-sec.  There  are  several  possible  chan- 
nels by  which  HF(v  = 3)  can  be  removed  by  H atoms;  unfortunately,  the  tra- 
jectory calculations*'  2 can  only  be  used  as  a rough  guide  since  the  London  - 
Eyring-Polanyi-Sato  (LEPS)  potential  energy  surfaces  used  differ  significantly 

7 

from  two  recent  quantum  mechanical  calculations  of  the  H-F-H  configuration. 

La  a previous  study,  ^ we  employed  the  laser-induced  fluorescence  tech- 
nique to  measure  the  room -temperature  removal  rate  of  HF(v  = 3)  and 
HF(v  * 2)  by  H atoms.  In  the  present  study,  the  measurements  were  extended 
to  200  K [for  HF(v  = 3)],  and  the  removal  rate  of  HF(v  = 3)  by  D was  measured 
as  a function  of  temperature.  With  a detailed  analysis  of  these  kinetic  data 
and  careful  consideration  of  the  F-H-H  thermochemistry,  the  following  con- 
clusions were  reached} 

1.  The  reaction  H + HF(v  = 3)  -•  H^  + F contributes  — 20% 

of  tiie  observed  removal  rate  at  295  K.  This  channel  must 
have  a small  positive  activation  energy. 


2. 


The  dominant  mechanism  for  HF(v  = 3)  removal  has 
a negative  temperature  dependence,  suggestive  of  com- 
plex formation,  attractive  interactions,  or  orbiting 
collisions.  The  possible  role  of  F-atom  exchange  is 
discussed. 


II.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 


The  basic  experimental  apparatus  has  been  previously  described.^  In 
the  present  experiments,  the  flow  tube  was  immersed  in  a stirred  low- 
temperature  slush.  A dry  ice -ethanol  mixture  was  used  for  the  200  K mea- 
surements, and  a water -ethanol  mixture  was  used  for  those  at  240  K.  Cold 
N2  was  slowly  passed  through  Cu  coils  at  the  bottom  of  the  bath  to  maintain 
the  slush.  The  temperature  of  the  slush  agreed  within  ~i  to  2°C  with  the 
reading  of  a thermocouple  mounted  in  the  center  of  the  flow  tube  downstream 
of  the  observation  point.  The  HF(v  = 3)  level  was  populated  by  sequential 
absorption  of  infrared  laser  photons  from  a multiband  pulsed  HF  laser.  This 

| 

level  was  monitored  by  3 -*  0 overtone  fluorescence  at  ~900  nm.  ^ Hydrogen 
(or  D)  atoms  were  created  by  a microwave  discharge  in  H2(or  D.,)  and  moni- 
tored by  isothermal  calorimetry.  ^ f 

Experiments  were  performed  at  total  pressures  of  ~3  Torr  with  partial  i 

pressures  of  ~1  x 10"3  Torr  HF,  ~3  X 10"2  Torr  H7,  and  <2  X 10'2  Torr 
H atoms;  He  made  up  the  balance.  1 


HI.  RESULTS 


r 


Within  the  precision  of  the  data,  the  3 -»  0 laser -induced  fluorescence 
traces  could  be  characterized  by  single  exponential  decays  for  the  first  one 
to  two  decay  times.  The  measured  fluorescence  traces  were  plotted  on  semi- 
log  paper,  and  the  exponential  decay  times  t were  determined.  In  each  ex- 
periment, the  decay  times  were  measured  with  the  microwave  discharge  on 
t and  with  the  discharge  off  t ,,  at  the  same  flow  rates.  The  measured 
decay  times  as  well  as  the  experimental  conditions  are  listed  in  Tables  I 
through  111.  The  overall  removal  rate  coefficient  of  HF(v  = 3)  is  designated 
k and  was  calculated  with  the  equation 

A(l/r)=klH]  (1) 

where  A(1/t)  = - t”^.  The  interpretation  of  this  k requires  discussion 

because  there  are  several  channels  for  HF(v  = 3)  removal. 

The  data  obtained  at  several  temperatures  are  plotted  in  Fig.  1.  ’ 

The  removal  rate  of  HF(v  = 3)  by  H is  faster  at  200  than  at  295  K.  Although  f 

the  temperature  range  is  small  and  the  scatter  is  somewhat  large,  the  data 

1 

imply  an  activation  energy  between  0 and  -750  cal/mol.  The  removal  rates  | 

I 

by  D atoms  are  ~10  to  20%  smaller  than  those  for  H atoms  at  both  295  and  1 

I 

200  K.  The  stated  uncertainty  in  the  deduced  rates  (30%)  reflects  both  system-  j 
atic  errors  in  the  isothermal  probe  and  flow-meter  calibrations  and  experi- 
mental scatter  in  the  data. 


T 

I 


Table  I . Removal  Rate  of  HF(v  = 3)  by  H Atoms  at  240  K 


Run 

^ total, 

[H]  X tO10. 

[H2]x  109. 

Ton' 

Toff’ 

A(I/t)  x 10-4. 
.1 

sec 

T. 

k X 10-13. 

t 

cm  Anol-sec 

Torr 

mol/ cm 

mol/cm3 

fxsec 

>10  ec 

*K 

1 

2.85 

9.  1 

2;0 

16.7 

136 

5.3 

240 

5.8 

2 

2.85 

9.7 

2.0 

12.  5 

136 

7.3 

238 

7.  5 

3 

2.85 

10.6 

2.0 

12.7 

136 

7.2 

236 

6.8 

4 

2.85 

6.94 

2.0 

18.4 

136 

4.7 

234 

6.8 

5 

2.85 

4.9 

2.0 

28.9 

136 

2.7 

235 

5.6 

6. 5 Average 


Table  n.  Removal  Rate  of  HF(v  = 3)  by  H Atoms  at  200  Ka 


Run 

Ptotal, 

Torr 

[H]X1010, 

mol/cm 

[H,]  x 109, 

c 3 

mol/cm 

T 

on, 

psec 

Toff, 

|isec 

A(1/t)x  10‘4, 
.1 

•cc 

T. 

*K 

. ,«-13 

k x 10  > 

cm^Anol- «ec 

1 

2.95 

7.2 

1.5 

13.8 

186 

6.7 

200 

9.3 

2 

2.95 

7.2 

1.5 

13.8 

186 

6.7 

200 

9.3 

3 

2.95 

6.2 

1.2 

17.  1 

186 

5.3 

200 

8.5 

4 

2.95 

6.4 

1.2 

17.6 

186 

5.  1 

202 

8.0 

5 

2.95 

8.8 

1.5 

11.2 

186 

8.4 

200 

9.6 

8.  9 Average 

1 

2.9 

7.4 

1.5 

11.8 

150 

7.8 

200 

11.0 

2 

2.9 

7.0 

1.  5 

11.4 

130 

8.0 

200 

11.0 

3 

2.9 

5.7 

1.2 

12.7 

130 

7.  1 

200 

13.0 

4 

2.9 

5.4 

1.2 

12.5 

130 

7.2 

200 

13.0 

5 

2.9 

6.4 

1.3 

13.7 

130 

6.5 

200 

10.0 

11.6  Average 

1 

2.95 

11 

2.2 

10.  1 

177 

9.3 

200 

8.6 

2 

2.95 

6.6 

1.5 

13.5 

177 

6.8 

200 

10.4 

3 

2.95 

6.7 

1.5 

14.5 

177 

6.3 

200 

9.5 

4 

2.95 

5.4 

1.5 

15.3 

177 

6.0 

200 

11.0 

9.9  Average 

*These  data  were  taken  on  three  separate  days. 


Table  III.  Removal  Rate  of  HF(v  = 3)  by  D Atoms  at  295  and  200  K 


Run 

Torr 

[DJ  xlO10. 

mol/cm^ 

[D,]  x 109, 

4 3 

mol /cm 

T , 
on 
l&iec 

To£f' 

)1MC 

A<1/t)X  10‘4. 
-1 

sec 

T. 

*K 

- 13 

kx  10  . 

cmVmol-iec 

1 

2.85 

8. 1 

3.2 

17.0 

146 

5.2 

295 

6.4 

2 

2.85 

9.0 

3.2 

14.5 

146 

6.2 

295 

6.9 

3 

2.85 

8.6 

3.2 

15.8 

146 

5.7 

295 

6.6 

4 

2.85 

9.5 

3.2 

15.8 

143 

5.6 

295 

5.9 

5 

2.90 

8. 1 

3.2 

18.5 

132 

4.7 

295 

5.8 

6 

2.90 

7.8 

3.2 

19.0 

132 

4.4 

295 

5.6 

7 

2.90 

6.0 

3.2 

27.0 

132 

3.0 

295 

5.0 

8 

2.90 

3.6 

3.2 

38.0 

132 

1.9 

295 

5.2 

5.9  Avenge 

9 

2.90 

8.5 

3.8 

13.7 

75 

6.0 

200 

7.0 

to 

2.9* 

6.0 

1.9 

14.8 

105 

5.8 

200 

9.7 

It 

2.90 

6.4 

t.9 

18.3 

105 

4.5 

200 

7.  1 

7. 9 Average 

a 

These  data  were  taken  on  two  separate  dcys. 
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The  disappearance  of  HF(v  = 3)  can  be  described  with  the  following 


kinetic  scheme. 


H + HF(v  = 3)  dL  H,  + F,  AH  = -617  cal/mol  (2) 

k 2 
-2 

k3 

H + HF(v  = 3)  -4.  H + HF(v  = 0,  1,  2),  AH  S -10,  356  cal/mol  (3) 

k4 

F + H2  -4-  HF(v  = 0.  1,  2)  + H,  AHS-9,  736  cal/mol  (4) 

k5 

HF(v  = 3)  + HF(v  = 0)  s*.  HF(v  = 2)  + HF(v  = 1),  AH  = 971  cal/mol  (5) 

k-5 

k6 

HF(v  = 3)  + HF  — HF(v  * 2)  + HF,  AH  = -10,  356  cal/mol  (6) 

k7 

HF(v  = 3)  + H2  — W HF(v  = 2)  + H2(v  = 0,  1)  (7) 

The  removal  of  HF(v  = 3)  can  occur  by  Reaction  (2),  which  in  the  backward 
direction  is  a channel  of  the  F + Hj  overall  reaction,  or  by  the  deactivation 
channel.  Reaction  (3).  The  latter  may  be  either  a nonreactive  inelastic 
collision  or  a reactive  collision  in  which  the  F atom  is  abstracted  by  the 


incoming  H atom.  With  no  H atoma  present,  the  decay  rate  of  HF(v  = 3) 
can  be  described  by 


r 


dfHFtv  = 3)1 
dt  [HF(v  = 3)] 


-<*5  * V [HF(v  . 0,]  ♦ k.5  tHF‘YH&^i)  • 1)1 


-k?  [H2]  - R 


(8) 


where  R represents  convection  and  diffusion  losses  out  of  the  observation 
volume  as  well  as  radiative  decay.  The  fluorescence  data  analyzed  in 
Tables  I through  III  decreased  monotonically,  indicating  that  the  right-hand- 
side  positive  (pumping)  term  was  always  smaller  than  the  negative  (removal) 
terms.  The  laser  could  be  detuned  or  the  [HF]  raised  so  that  this  was  not 
true,  i.  e.  , the  fluorescence  increased  initially  and  then  decayed.  With  the 
deactivation  rates  listed  in  Table  IV  and  typical  concentrations  of 
HF~5  X 10"11  and  H2~l.  5 X 10'9  mol/cm3,  (kg  + kfe)  [HF]  ~2  X 103  sec"1 

2-1  3 

and  k^  [H2]  **3  X 10  sec  These  rates  and  an  estimate  of  R ~ 2 x 10 

-1  3-1 

sec  indicate  a total  decay  rate  of  ~4  X 10  sec  and  a decay  time  of 

'"’250  pace,  which  are  in  reasonable  agreement  with  the  observed  decay  times 

at  room  temperature. 

When  the  laser  was  tuned  to  maximize  the  HF(v  = 3)  fluorescence,  the 
signal  decayed  exponentially  for  one  to  two  decay  times  when  no  H atoms 
were  present.  At  longer  times,  the  decay  rate  decreased,  probably  because 
the  pumping  term  k_g  becomes  increasingly  important  at  the  smaller 
HF(v  a 3)  concentrations.  Hydrogen1*  does  not  contribute  greatly  to  the 


Table  IV.  Reaction  Rate  Coefficients 


Rate 

Coefficient 

Assumed  Value, 
cm'Vnol- sec 

Reference 

k-2 

4.2  X 1012 

10,  see  text 

k2 

13 

1.37  x 101J 

See  text 

k4  +k-2 

13 

1 . 5 x 101J 

10 

k5  + k6 

13 

2.9  X 10 

11 

k7 

2.0  x 1011 

12 

decay  rate  even  though  the  [H2J  * 30  [HF];  the  back  reaction.  Reaction  (-7), 
is  insignificant  compared  with  Reaction  (-5). 

When  H atoms  are  present,  two  additional  terms  must  be  added  to 
the  right-hand  side  of  Equation  (8).  They  are 


[F]  [H2] 

(k2  + k3)  ^ + k-2  [HF(v  = 3)] 


The  sum  of  these  two  terms  has  a large  negative  value,  compared  to  the  other 
terms  in  Equation  (8),  since  the  HF(v  = 3) decayed  typically  five  to  ten  times 
faster  with  the  H atoms  than  without.  Therefore,  the  disappearance  rate  of 
HF(v  = 3)  can  be  approximated  by 

- (k5  + k6)  [HF(v  = 0)]  - k?  [H2]  - R (10) 


If  the  second  term  on  the  right-hand  side  of  Equation  (10)  is  small,  the  decay 
is  essentially  exponential.  The  change  in  decay  rate  produced  by  the  H atoms 
can  be  described  by 


A (t)  ' K + k3  - T k7)  [“1  * <k2  ♦ k3»  im 


This  change  in  decay  rate,  if  evaluated  at  time  equal  to  aero,  is  always 
given  by  Equation  (11).  However,  if  k2  » kj,  the  effect  of  the  pumping  term 


9 


i 


in  Equation  (10)  must  be  taken  into  account  if  the  decay  rates  are  evaluated 
at  longer  times. 

It  is  shown  in  the  Appendix  that  the  preceding  equations  are  equivalent 
to  those  describing  the  V-V  and  V-R,  T relaxation  of  an  initially  excited 
molecular  species  by  a chaperone  molecular  species.  The  concentration  of 
HF(v  = 3)  can  be  described  by  the  sum  of  two  exponentially  decaying  terms. 

[HF(v=  3)]  = A exp(-Xjt)  + B exp(-X2t)  (12) 

The  solution  is  given  in  the  Appendix  to  a somewhat  simplified  set  of 
reactions.  The  values  of  the  two  inverse  decay  times  X^  and  normalized 
with  the  quantity  Q + N are  plotted  in  Fig.  2.  They  were  calculated  for 
lH2]  = 3 x IQ'9  mol/cm^  with  the  rate  coefficients  of  Table  IV  and  three 
values  of  k^.  In  summary,  Q = k2[H],  and  N is  the  sum  of  all  other 
deactivation  (negative)  terms  in  Equation  (10).  Q + N is  always  the  initial 
deactivation  rate  of  HF(v  = 3);  the  initial  deactivation  rate  should  not  be 
confused  with  Xj.  F atoms  are  produced  by  Reaction  (2)  until  an  equilibrium 
is  established  with  Reaction  (4),  which  removes  them.  This  equilibrium  is 
established  quickly  when  H2  is  large  compared  to  H.  Thereafter,  HF(v  = 3) 
decays  with  the  rate  of 

x2-(s77tt;)q  + n (,j> 


8+V  N+D  N+D 


► 


where  k^/(k^  +k  ^)  hae  been  determined^  to  be  ~0.  7.  This  decay  rate  is 
at  most  only  30%  slower  than  the  initial  decay  rate  and  more  nearly  equal  to 
it  if  any  substantial  deactivation  occurs  by  Reaction  (3),  i.  e.  , if  5 k^. 

At  larger  H atom  concentrations,  [H]  ~ [H^],  a larger  fraction  of  the 
HF(v  = 3)  has  to  dissociate  before  the  pumping  term  in  Equation  (10), 
k_2{[F]  [H2J/[HF(v  = 3)]},  becomes  significant  compared  to  the  deactivation 
terms,  (k^  + k^)  [H].  Therefore,  the  initial  decay  rate,  Q + N,  should  per- 
sist for  a larger  proportion  of  the  HF(v  = 3)  decay.  The  two  decay  rates 
should  be  the  most  apparent  and  the  most  easily  resolved  when  A~  B in 
Equation  (12)  and  when  the  two  decay  rates  are  quite  different.  For 
A/(A  + B)  = 0.  5,  X2  is  only  a factor  of  3.  3 slower  than  Q + N,  the  initial 
decay  rate,  if  k2  9P  k^  (see  the  curve  for  N = 0,  Fig.  2).  If  k^  i k2>  there 
is  even  less  of  a spread  between  the  two  decay  rates.  The  difficulty  of 
determining  whether  a measured  trace  decays  as  a single  exponential  or  as 
the  sum  of  two  exponentials  is  demonstrated  in  Fig.  3,  where  a calculated 
fluorescence  trace  described  by  Equation  (12)  is  plotted  with  the  values  for 
\j,  *2'  an<*  A/(A  + B)  taken  from  Fig.  2.  The  circles  represent  the  same 

trace  except  that  0.  05  X I was  subtracted  to  demonstrate  the  effect  of  the 

o 

typical  uncertainty  in  the  baseline  of  such  experimental  data.  The  result 
is  that  the  circles  can  be  fitted  within  2%  of  the  full-scale  reading  with  a 
single  exponential  decay  rate.  Note  that  the  approximate  fit  to  the  circles 
has  a decay  rate  within  10%  of  the  initial  decay  rate  of  the  theoretical  trace. 
The  precision  of  the  data  in  these  experiments  is  not  sufficient  for  two 
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decay  rates  to  be  resolved,  much  less  quantitatively  established.  The  two 
decay  rates  would  be  much  more  distinguishable  if  a larger  fraction  of  the 
F + reaction  were  produced  in  HF(v  = 3)  than  the  currently  estimated 
fraction  of  0.  28,  i.  e.  , if  Reaction  (4)  were  slower  relative  to  Reaction  (-2). 

In  the  light  of  these  uncertainties,  we  can  only  interpret  the  measured 
values  of  k in  Equation  (1)  as 


k = A k^  + k^ 


(14) 


with  A between  0.  7 and  1. 

The  rate  coefficient  for  Reaction  (2)  has  not  been  measured.  However, 
it  can  be  estimated  on  the  basis  of  the  overall  reaction  rate  for  F + H^,  the 
relative  vibrational  distributions,  and  the  equilibrium  constant  - 

k^/k  Cohen  and  Bott*^  reviewed  the  data  for  the  F + reaction  and 
recommend  a rate  coefficient  of 


1 


< 


kF+H  = 2.  3 X 1014exp(-1600/RT) 


(15) 


13  3 

with  a value  of  (1.  5 ± 0.  5)  X 10  cm  /mol-sec  at  295  K.  In  the  same 

review,  Cohen  and  Bott  recommend  a value  of  0.  28  for  the  branching  fraction 

into  the  v = 3 level,  i.  e.  , k_£/(k_2  + k^)  = 0.  28,  or  a value  of  4.  2 x 10*2 

cm  /mole-sec  for  k^.  A value  of  ^ = 3.2  at  295  K can  be  calculated 

13 

from  the  JANAF  thermodynamic  data  (Table  V).  The  uncertainty  in  this 
value  stems  largely  from  the  uncertainty  in  the  bond  dissociation  energy 
(BDE)  of  HF  since  the  exothermicity  of  Reaction  (2)  depends  directly  on  the 
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Table  V.  JANAF  Thermodynamic  Data 
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Species 

AHi-,98 

kcal/mol 

H 

52.  100 

D 

53.252 

HF(v  = 0) 

-65.  140 

HF(v  = 3) 

-32.623 

Hz 

0.000 

HD 

0.076 

F 

18.860 

I 


c 


<•! 


I 

j 

\ 

\ 

bond  dissociation  energies  of  HF  and  and  the  spectroscopically  determined 
energy  levels  of  HF.  The  exothermicity  of  Reaction  (2)  at  0 K is  simply 

AE  = BDE(H2)  - BDE(HF)  + EHF(3)  " EHF(0)'  where  EHF(3)  " EHF(0)  i-  the 

energy  difference  between  HF(v  = 3)  and  HF(v  = 0).  The  JANAF  thermo- 

14 

dynamic  data  (Table  V)  are  based  upon  a value  of  135, 120  ± 300  cal/mol 

for  BDE(HF)  at  0 K.  A more  recent  measurement  by  Di  Lonardo  and 
1 5 

Douglas  gave  a value  of  135,  274  ±170  cal/mol,  which  agrees  with  the 

earlier  measurement.  An  uncertainty  of  200  cal/mol  produces  an  uncertainty 

of  a factor  of  1. 4 in  the  calculated  value  of  K2  2 295  K [a  larger  value 

for  BDE(HF)  gives  a smaller  value  for  K2  _2], 

With  the  value  of  3.  2 for  K2  _2  and  the  recommended  value  for  k_2, 

13  3 

a rate  coefficient  of  k2  = 1.  37  x 10  cm  /mol-sec  is  obtained,  compared  to 
13  3 

(6.  3 ± 1.  5)  x 10  cm  /mol-sec  obtained  in  the  experiments  for  the  total 

removal  rate  at  295  K.  ^ These  two  values  differ  by  a factor  of  4.  6.  The 

uncertainty  of  a factor  of  1.  5 in  the  overall  F + H2  rate*^  and  the  fraction 

of  the  reaction  going  into  v = 3 combined  with  the  factor  of  1.  4 uncertainty 

in  K2  _2  gives  a total  uncertainty  of  a factor  of  2 in  the  estimated  value  of 

1 3 

k2<  The  estimated  upper  limit  of  2.  7 x 10  for  k2  is  a factor  of  1.  8 lower 

13  3 

than  our  estimate  of  the  lower  limit  of  4.  8 x 10  cm  /mol-sec  for  the 
removal  rate. 

Although  there  have  been  no  prior  experimental  studies  of  Reaction  (2), 

Wilkins2  has  performed  trajectory  calculations  on  an  LEPS  potential  energy 

12  3 

surface.  He  obtained  a value  of  4,  2 X 10  cm  /mol-sec  for  k2,  compared 


-2 


’i'.vnl 


lo  the  value  of  1.  37  x 10  , which  was  calculated  from  the  equilibrium 

constant  and  the  value  of  k_g  recommended  by  Cohen  and  Bott.  Xn  a sepa- 
rate calculation,  ^ k g = 2.  66  x 10*^  cm3/mol-sec  was  obtained;  the  recom 
mended  value  is  4.  2 x 10  * cm  /mol-sec.  Wilkins'  two  calculated  rates 
yield  a value  of  1.  6 for  K g _£>  a factor  of  2.  0 lower  than  that  calculated 
from  the  JANAF  data.  However,  a factor  of  2 is  probably  within  the  accu- 
racy of  the  trajectory  calculations.  Wilkins'  calculations  suggest  that  kg 
may  be  <10%  of  the  total  measured  removal  rate  of  HF(v  = 3). 
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The  temperature  dependence  of  Reaction  (2)  can  be  estimated  from 

1 7 

existing  data.  Coomb e and  Pimentel  measured  the  temperature  dependence 
of  the  fraction  of  the  F + Hg  reaction  going  into  v = 3,  relative  to  the  fraction 
yielding  v = 2.  They  found  that 

k-2 

= 0. 39  exp{  1 17/RT)  (16) 

Approximately  55%  of  the  F + Hg  reactive  collisions  yield  HF(v  = 2);  there- 
fore, the  temperature  dependence  of  k^(v  = 2)  can  be  approximated  with  that 
of  the  overall  reaction.  Combining  Equations  (15)  and  (16)  produces  the 
expression 

k « = 4.  9 X 1013  exp  (-14Q3/RT)  (17) 


u 


The  equilibrium  constant  determined  from  the  JANAF  data  can  be  described 
between  200  and  400  K as  follows: 

K2  _2  = 0.  338  T°*  179  exp(723/RT)  (18) 

Thus,  an  expression  for  k2  can  be  derived  from  Equations  (17)  and  (18) 


k2  = 1.  7 X 1013  T°*  179  exp(-760/RT) 


(19) 


Written  as  a simple  Arrhenius  expression,  k2  has  an  effective  activation 

energy  of  850  cal/mol,  which  is  unlike  the  negative  activation  energy  of  the 

measured  removal  rate.  The  calculated  value  of  k2  at  T = 200  K is 
13  3 

~0.  7x10  cm  /mol-sec,  roughly  7%  of  the  measured  removal  rate  of 

13  3 6 

10  X 10  cm  /mol-sec.  As  earlier  work  concluded,  the  reaction  to  form 

H2  + F is  insufficient  to  explain  the  removal  of  HF(v  = 3)  by  H. 

C.  REMOVAL  OF  HF(v  = 3)  BY  P ATOMS 

13 

The  measured  removal  rates  of  HF(v  = 3)  by  D atoms  (5.9  x 10 
cm3 /mol-sec  at  295  K and  7.  9 x 10*  3 at  200  K)  were  within  20%  of  the 
slightly  faster  rates  measured  for  HF(v  = 3)  removal  by  H atoms.  The 
removal  of  HF(v  ■ 3)  by  D atoms  can  be  described  with  the  following  set  of 
reactions: 


HD  + F,  AH  = -1691  cal/mol 


D + HF(v  = 3) 


20 


(20) 


D + HF(v  = 3)  -iift  D + HF(v  = 0,  1,  2),  AH  < -10,  356  cal/mol  (21a) 


H + DF 


(21b) 


and  possibly 


F + Dz— *D  + DF 


For  the  conditions  of  the  present  experiment,  the  concentration  of  HD  was 
always  small  compared  to  that  of  so  that  Reaction  (-20)  can  be  neglected, 
and  the  HF(v  = 3)  can  be  expected  to  decay  exponentially.  The  difference  in 
the  decay  rates  with  and  without  the  microwave  discharge  is  simply 


4 (?)  * <k20  * k2l>  M 


where  the  term  -1/2  was  neglected  with  respect  to  (k2Q  + k21). 

2 

The  rate  of  HF(v  = 3)  relaxation  by  D^,  _ D , was  found  in  a separate 

18  12  3 ^ 

study  to  have  a value  of  1 x 10  cm  /mol- sec  at  295  K. 

19 

Berry  determined  the  overall  rate  for  HF  formation  by  the  F + HD 
reaction  to  be  slower  than  that  for  F + H2  by  a factor  of  2.  5 at  room  tem- 
perature. He  also  determined  that  Reaction  (-20)  accounted  for  10%  of  the 
overall  F + HD  -•  HF  + D reaction.  With  Berry's  results  and  the  value  of 
1.  5 x 10*3  cm  /mol-sec  that  Cohen  and  Bott*®  recommended  for  the  F + H2 
rate  coefficient  at  299  K, 


k 2q  = 6.  0 x 10**  cm^/mol-sec  at  295  K 


can  be  calculated.  The  equilibrium  constant  I^q  _2q  is  35  at  295  K 

according  to  a calculation  based  on  the  JANAF  thermodynamic  data 

13  3 

(Table  V).  Therefore,  k^Q  has  a value  of  2.  1 x 10  cm  /mol-sec  at  295  K, 

1 3 

somewhat  larger  than  the  value  of  k^  = 1.  37  X 10  for  the  comparable 

reaction  for  H instead  of  D.  This  value  of  k^Q  is  ~2.  8 smaller  than  the 

13  3 

total  removal  rate  of  5.  9 X 10  cm  /mol-sec  measured  in  the  experiment. 

The  value  of  k^Q  is  uncertain  in  the  same  proportion  as  k ^ since  they 
are  both  derived  from  the  value  of  k It  has  an  additional  uncertainty 
associated  with  Berry's  measurements  of  the  F + HD  — HF  + D reaction 
rate  relative  to  that  for  F + — HF  + H and  the  fraction  going  into  HF(v  = 3). 

He  obtained  distribution  numbers  for  the  pumping  reactions  only  slightly 
different  from  the  accepted  values.  ***  There  have  been  no  measurements  of 
the  temperature  dependence  of  Reaction  (20);  therefore, there  is  no  guide  to 
the  rate  at  200  K.  It  should  be  mentioned  that  Wilkins^  made  trajectory 
calculations  on  an  LEPS  potential  energy  surface  for  Reaction  (-20)  and 
obtained  a rate  coefficient  of 


k ,q  = 1.35  x 10l  5 exp  (-1628/RT)  cm  /mol-sec 


which  has  a value  of  8.  4 x 10**  cm  /mol-sec  at  295  K,  compared  to  the 
113 

value  of  6. 0 X 10  cm  /mol-sec  calculated  from  Berry's  data  and  the  recom- 
mended rate  for  F + H^.  Equation  (24)  has  an  activation  energy  very  close 


to  the  exothermicity  of  Reaction  (20).  Therefore,  it  ie  inferred  from 
Wilkins'  rate  for  k and  the  JANAF  thermodynamic  data  that  the  activation 
energy  of  k^Q  is  essentially  zero. 

D.  COMPARISONS  WITH  OTHER  DATA 

The  only  previous  experimental  study  of  HF(3)  removal  by  H atoms 

3 

was  a flow-tube  study  by  Kwok  and  Wilkins,  However,  in  their  experiments, 
there  were  equal  parts  of  HF  and  H since  the  HF(v)  was  produced  by  reacting 
F atoms  with  H^.  The  decay  rates  were  measured  at  various  H^  concen- 
trations and  then  extrapolated  back  to  [H^]  =0,  where  the  decay  rate  should 
be  the  result  of  processes  involving  HF  and  H and  spontaneous  emission. 
Their  measurements  extrapolated  to  [H2]  = 0 yielded  1.8  X 10  cm  /mol- sec 
for  the  rate  of  removal  of  HF(3)  by  H and  HF.  This  value  is  a factor  of  3.  5 
slower  than  the  present  results,  even  though  it  contains  the  additional  con- 
tribution of  HF  self- relaxation.  On  the  other  hand,  Kwok  and  Wilkins 
deduced  a much  faster  rate  coefficient  (by  a factor  of  20)  for  the  deactivation 
of  HF(v  a 2)  by  H atoms  than  our  measurements  indicate.  ^ Since  the  flow- 
tube  values  are  slower  for  v * 3 and  faster  for  v = 2,  it  is  possible  that  the 
effects  of  HF-HF  V-V  coupling  have  to  be  taken  into  account.  The  inter- 
pretation of  those  flow-tube  experiments  is  still  under  investigation. 

Two  theoretical  studies*’  * of  the  relaxation  of  HF  by  H atoms  in 
Reaction  (3)  have  been  made;  both  were  Monte  Carlo  classical  trajectory 
calculations.  Thompson's1  calculations  were  performed  at  temperatures 
a 600  K,  but  his  results  indicate  that  the  deactivation  rate  increased 


*4 


2 

approximately  as  the  vibrational  level  v.  Wilkins  calculated  the 

deactivation  rates  of  several  vibrational  levels  of  HF  by  H atoms  with  a 

semi-empirical  LEPS  potential  energy  surface.  The  initial  calculations  were 

performed  with  a 1500  cal/mol  barrier  height  for  H-F-H,  but  his  results 

3-5 

for  the  v = 1 level  were  faster  than  experimental  results;  therefore,  the 

effect  of  barrier  height  was  examined.  ^ Calculations  with  barriers  of 

12 

1500,  2500,  and  3500  cal/mol  gave  deactivation  rates  of  2.  5 x 10  , 

3. 1 X 10**,  and  3.  1 x 10***  cm^/mol-sec  for  the  first  vibrational  level  at 
room  temperature;  the  value  of  3.  1 X 10**  agreed  with  the  experimental 
data.  In  each  case,  the  deactivation  rate  of  HF(v  = 3)  was  ~ 6 times  faster 
than  that  for  v = 1,  and  the  final  vibrational  states  were  approximately 
equally  distributed  over  v = 0,  1,  and  2,  with  about  equal  proportions  of 
reactive  and  nonreactive  deactivation. 

7 

Recently,  an  ab  initio  calculation  by  Bender,  Garrison,  and  Schaefer 

indicated  that  the  H-F-H  surface  has  a barrier  height  of  ~40  kcal/mol. 

8 9 

Preliminary  results  of  a similar  calculation  by  Wadt  and  Winter  ' also 

indicate  a high  barrier  of  36  ± 4 kcal/mol  and  differ  markedly  from  an 

LEPS  angular  dependence.  Wilkins^*  performed  a trajectory  calculation 

on  an  LEPS  surface  with  a 40  kcal/mol  barrier  and  obtained  a deactivation 
12  3 

rate  of  2.  5 x 10  cm  /mol-sec.  Thus,  all  of  the  trajectory  calculations 
gave  significantly  smaller  values  for  k^  than  the  values  of  (k^  + kj)  obtained 
in  the  present  study.  The  calculations  also  exhibit  positive  activation 
energies,  compared  to  the  negative  activation  energy  of  the  present 
measurements. 


The  vibrational  level  v = 3 of  HF  is  38.  3 kcal/mol  above  the  bottom 


of  the  well  on  the  potential  energy  surface,  which  is  comparable  to  the 

7 

barriers  to  F-atom  transfer  estimated  by  Bender  et  al.  and  by  Wadt  and 
Winter.  ^ Smith  and  Wood^  investigated  the  relaxation  of  vibrationally 
excited  molecules  when  an  atom  exchange  is  possible.  They  found  that  vibra- 
tional excitation  above  the  barrier  not  only  permits  atom  exchange  with  con- 
sequent loss  of  vibrational  energy  but  also  allows  multiple  barrier  crossings, 
which  increase  the  possibility  for  conversion  of  vibrational  energy  to  rotational 
or  translational  energy,  regardless  of  whether  or  not  reaction  occurs. 


A 
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V.  CONCLUSIONS 


* 


The  experimental  results  are  not  precise  enough  to  make  a complete 
allocation  of  the  measured  removal  rate  to  the  possible  removal  channels. 

If  the  removal  rate  were  the  result  of  Reaction  (2),  die  fluorescence  traces 
would  exhibit  the  sum  of  two  exponential  decay  terms.  These  decay  terms 
are  related  in  such  a way  that  extremely  precise  data  would  be  required  to 
extract  the  separate  decay  rates.  The  present  data  fit,  and  are  analyzed  in 
terms  of  a single  exponential  decay  rate  k = Ak2  + k^,  where  0.7  < A<  1. 

Reaction  (2)  contributes  substantially  to  the  measured  decay  rate. 
Accepted  thermodynamic  and  kinetic  data  for  the  H + HF  system  were  used 
to  estimate  k2/k  = 0.2  \ at  295  K.  Similarly.  k2(J/(k20  +k2J)  = 0.4 

has  been  estimated  for  the  D + HF(v  = 3)  deactivation  at  295  K.  The 
largest  uncertainty  in  the  thermodynamic  data  for  Reaction  (2)  is  in  the  bond 
dissociation  energy  of  HF.  An  uncertainty  of  200  cal/mol  in  the  value  pro- 
duces an  uncertainty  of  a factor  of  1.4  in  the  calculated  equilibrium  constant 
for  the  reaction  at  295  K. 

The  removal  rate  for  HF(v  = 3)  by  H atoms  is  ~ 100  times  faster  than 
that  for  HF(v  = 2).  Unless  Reaction  (2)  represents  >95%  of  this  rate,  which 
is  highly  improbable,  the  V -*  R,T  deactivation  of  HF{v  = 3)  must  be  much 
faster  than  the  V -*  R,  T deactivation  of  HF(v  = 2).  A collision  of  HF(v  = 3) 
with  an  H atom  can  easily  have  sufficient  energy  to  result  in  the  formation 
of  H2>  Md,  if  recent  quantum  mechanical  calculations  are  correct,  there  may 
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be  sufficient  energy  for  the  incoming  H atom  to  abstract  the  F atom.  Even  if 

22 

these  atom  transfers  do  not  occur.  Smith  and  Wood  found  greatly  enhanced 

vibrational  deactivation  probabilities  in  theoretical  studies  of  similar 

triatomic  systems  in  which  atom  exchanges  were  energetically  possible. 

The  T~*  temperature  dependence  observed  for  the  rate  coefficient  is 

23  24 

very  similar  to  that  observed  for  the  relaxation  of  HF(v  = 1)  ' and 

HCl(v  = 1).  25  Zittel  and  Moore25  discussed  the  HC1  and  HBr  results  and 
their  possible  explanation  in  terms  of  attractive  potentials,  multiple  encounters, 
long-lived  bimolecular  collisions,  and  complex  formation.  These  mechanisms 
suggest  increasing  vibrational  relaxation  rates  at  decreasing  temperature 
and  may  apply  to  the  deactivation  of  HF(v  = 3)  by  H atoms.  The  present 
results  do  not  agree  well  with  published  trajectory  calculations,  and  additional 
calculations  made  with  approximations  to  a more  realistic  surface  would  be 
very  useful. 
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APPENDIX 


F + H2  -i  HF(0,  1,2)  + H (A- 3) 


The  differential  equations  that  describe  these  reactions  under  the  conditions 
of  the  present  experiments  are  equivalent  to  those  for  the  vibrational  energy 
relaxation  processes: 


4 


Y*  = [F]  = Cfexpf-Kjt)  - exp(-\2t)]  (A-8) 

23 

The  two  decay  times  Xj  and  can  be  expressed  as 

X1  + X2=Q+P  + N + K CA-9) 

Xj  X X2  = QK  + PN  + NK  (A-  1C) 

where  Q = k ^ [H],  P = k_2fH2J,  N = [H],  and  K = k^  [H2J. 

From  Reactions  A- 1 and  A-2,  it  is  seen  that 

dfHFgfl  = - [HF(3)J  (Q  + N)  + [F]P  (A-  1 1) 

dt 

Equations  (A- 7)  and  (A-8)  can  be  substituted  into  Equation  (A- 11).  After 
collecting  the  terms  containing  exp(-Xjt)  and  exp(-X2t),  two  equations  are 
obtained  which  can  be  solved  for  A/ (A  + B)  to  yield: 

Q+N-X. 

A/(A+B)=— r r— £ (A- 12) 

*1  ' *2 

The  initial  exponential  decay  rate  is  always  Q + N since  [F]  is  initially 
sero  in  Eq.  (A-ll).  Therefore,  it  is  convenient  to  normalize  Xj  and  X2 
with  Q + N.  Calculations  have  been  performed  with  the  rate  coefficients 
listed  in  Table  IV  for  the  conditions  of  [H2J  = 3.  0 x 10"9  mol/cm3  and 
various  H atom  concentrations.  The  normalised  values  for  X j,  X2,  and 
A/(A  + B)  have  been  plotted  versus  [h]/[H21  in  Fig.  2 for  three  values  of 
N * k3  [H] . 
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